Quantum kinetic theory for fermions in temporally 

varying backrounds 



Matti Herranen, Kimmo Kainulainen and Pyry Matti Rahkila 

University of Jyvdskyld, Department of Physics, P. 0. Box 35 (YFL), 
FIN-40014 University of Jyvdskyld, Finland 

O . and 

O 

Helsinki Institute of Physics, P.O. Box 64, FIN-OOOI4 University of Helsinki, Finland. 



'~~^- Abstract 

^ ■ 

jiLi' We derive quantum kinetic equations for fermions in a homogeneous time-dependent back- 

(-1 ■ ground in presence of decohering colhsions, by use of the Schwinger-Keldysh CTP-formahsm. 

The quantum coherence (between particles and antiparticles) is found to arise from new 

spectral solutions for the dynamical 2-point correlation function in the mean field limit. The 
if^ • physical density matrix p and its dynamics is shown to be necessarily dependent on the ex- 

^ ■ trenous information on the system, and expressions that relate p to fundamental coherence 

functions and fermionic particle and antiparticle numbers are derived. For an interacting 
t^^ ■ system we demonstrate how smooth decoherence effects are induced by collisions. As special 

QQ ■ applications we study the production of unstable particles during the preheating stage of 

O ' the inflation and an evolution of an initially quantum p towards a statistical limit including 

decoherence and thermalisation. 
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1 Introduction 

We study the quantum dynamics of fermions in a homogeneous but temporally varying 
background field including collisions with a thermal background. These conditions are ap- 
propriate to model for example the particle production during preheating at the end of 
infiation or during cosmological phase tranistions P,[2], as well as baryogenesis during pre- 
heating [3], or coherent baryogenesis [Ij. They are relevant also for neutrino oscillations in 
the early universe [5], or just for generic studies of thermalisation of quantum systems [6]. 
The formalism we will be developing here can also be modified for a treatment of static 
problems with planar symmetry [TJ E] . In this form the resulting kinetic equations will be of 
interest for problems involving quantum refiection such as electroweak baryogenesis [HI ^U\ . 
For related studies of quantum transport equations for electroweak baryogenesis starting 
from similar theoretical foundation see e.g. [TT1[T2]. 

By use of the CTP-formalism [131 [H] we will set up quantum kinetic equations for 
the fermionic 2-point function G"^. We find that these equations admit a rich structure of 
spectral solutions including the expected mass-shell solutions for particle and antiparticle 
excitations, but also a new class of solutions living in k^ = 0-shell in the phase space [7j. We 
interpret the arbitrary weight functions on these shells as describing the out-of-equilibrium 
particle and antiparticle numbers (mass shells) and the quantum coherence between particles 
and antiparticles {ko = 0-shell). New coherence solutions are shown to be excluded from 
the spectral function by the spectral sum-rule, which indicates that they are not part of the 
kinematic phase space, although they necessarily occur in the dynamical function G'^. These 
solutions are also eliminated from G^ by the KMS-conditions in the thermal limit. As a 
consequence of the singular shell structure, an integration procedure is needed in order to 
define a physical density matrix in terms of the original 2-point function. We show how this 
procedure necessarily involves specifying precisely the amount of information on the system, 
and derive an evolution equation for a density matrix relevant for a spatially homogenous 
case, including quantum coherence. We then introduce the interactions and show that the 
spectral structure for the phase space including the coherence shell solutions survives in the 
quasiparticle limit with the interactions. We set up dynamical equations for the physical 
density matrix including the interactions and compute explicit expressions for the collision 
integrals in the case of a simple model Lagrangian describing decays and inverse decays. 
The interaction terms are shown to contain the usual collision terms that push the mass- 
shell functions towards thermal limit, but also other collision terms that tend to bring the 
quantum coherence functions to zero. The loss of coherence does not happen instantly, like a 
collapse of a wave function, but smoothly over a characteristic time scale set by the strength 
of the interactions at the shell ko = 0, in close analogy to the damping of coherence in the 
case of neutrino oscillations [^ . 



In our approach we define the fermionic number density in the same way as in thermal 
field theory, that is, as an dimensionless real- valued function living on the positive mass-shell 
of the spectral form 2-point correlator. With this definition we see that in thermal limit the 
number density is indeed the standard Fermi-Dirac distribution. We will show that our 
definition for the particle number agrees with that of ref. [2], where it was derived using the 
operator methods and Bogolyubov transformation to diagonalize the fermionic Hamiltonian. 
We also derive expressions for the energy density and the pressure. The latter is shown to 
differ from the statistical pressure, but the statistical pressure is retrieved for any realistic 
measurement that averages out the quantum oscillations. As applications of our formalism 
we calculate particle number production during fermionic preheating including finite decay 
width for the heavy particles produced. We will see that the decoherence induced by the 
decays can have dramatic effect on the particle number evolution. We also show explicitly 
how an initially highly correlated out-of-equilibrium density matrix relaxes to a thermal 
equilibrium as a result of collisions. 

This paper is organized as follows: In section [2] we give a brief intoduction to Shwinger- 
Keldysh CTP-formalism in order to write down the fermionic Kadanoff-Baym (KB) transport 
equations for 2-point correlation functions. In section [3] we study the noninteracting KB- (or 
Dirac) equation and find out the nontrivial phase-space shell structure including the new 
ko = 0-shell. In section H] we define a nonsingular physical density matrix, as a weighted 
integral over the initial singular 2-point function and show that its evolution is heavily 
dependent on the extrenous information we have on the system. The material in the sections 
|3l|4]partly overlaps with the derivation in the companion paper [7], but we include a shortened 
discussion here for completeness. In section [5] we will compute the particle number density, 
energy density and pressure in terms of the density matrix. In section E] we generalize our 
kinetic equations to include collisions. We study the particle production at the preheating 
and the decoherence phenomenon in sections 17.11 and 17.21 and finally section [H] contains our 
conclusions. 

2 General fermionic CTP-formalism 

The basic object of interest in this paper is the fermionic 2-point Wightmann function defined 



as: 



iG<f^{u,v) = {ii)p{v)i)a{u)) = Tr {p i)p{v)i)^{u)] 



where p is some unknown quantum density operator describing the properties of the system. 
In a non-interacting theory G^ decouples from other n-point functions and the dynamical 
equation it satisfies is equivalent to the ordinary Dirac equation. This problem was studied 
carefully in ref. [7j. However, here we wish to include interactions, and so it is necessary to 
work in the framework of the quantum field theory. The QFT formalism that is well suited 
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Figure 1: Schwinger-Keldysh path in complex time. 

for the study of "in-in" -correlator^ like (^ in possibly out-of-equilibrium conditions is called 
Schwinger-Keldysh or Closed Time Path (CTP) formalism [13l[ll]. In that formalism one 
defines a path ordered 2-point function on a complex Keldysh time-path (Dirac indices are 
suppressed) : 

tGc{u,v) = {Tc[iP{u)^p{v)]), (2) 

where Tc defines time ordering along the Keldysh contour C, which starts at some to, often 
taken to be at — cxd, goes to +cx3, and then back to to (see Fig. [1]). It can be shown for 
example by using the two-particle-irreducible (2PI) effective action techniques P^ [T3j that 
the 2-point function Gc{x,y) obeys the contour Schwinger- Dyson equation: 

Gc{u,v) = G'c{%v)+ fd^z,fd''z2G''c{u,z,)Ec{zuZ2)Gc{z2,v), (3) 

Equation ([3]) is formally expressed in Fig. [21 where the thin lines correspond to the free 
particle (tree level) propagator G°, and the thick lines to the full propagator Gc- The filled 
ellipsis represents the self-energy function Sc, whose precise form depends on the model 
Lagrangian and a truncation scheme. Ec couples Gc to an infinite (BBGKY-) hierarchy of 
equations for higher (up to infinite) order Green's functions, and some approximation scheme 
is needed to truncate this hierarchy to obtain the closure. In the weak coupling limit it will 
be natural to do this by substituting all higher than 2-point functions by their perturbative 
expressions. However, we can learn a lot about the structure of the SD-equations ([3]) without 
ever making any reference to the explicit form of S. Multiplying Eq. (^ by the inverse of 
the free particle propagator (G^)"^ and integrating over the connecting variable zi one finds 



f d'zG'^ciu, z)-'Gc{z, v) = 6c{u-v)+ I d'zJ^du, z)Gc 



[z,v), 



(4) 



where 5c{u — v) = Sc{uq — Vq)6^{u — iT) is a contour time delta-function. The complex time 
Green's function in (I21I11) can be conveniently decomposed in four different 2-point functions 



^With "in- in" -correlator we mean that the matrix elements are expectation values {in\A\in), in contrast 
to the traditional QFT transition amplitudes {out\A\in) with different in and out states 
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Figure 2: A generic form of a Schwinger-Keldysh equation for the 2-point function G^. 

with respect to usual real time variable: 

tG<{u,v) = -iG+~{u,v) = {^{v)i/j{u)) 

iG>{u,v)= iG-+{u,v) = {ij{u)tlj{v)) 

iGpiu^v)^ iG^~^{u,v) = 6{uq — vq)G^{u,v) — 6{vq — uq)G^{u,v) 

iGp{u,v) = iG {u,v) = 9{vo — uq)G^{u,v) — 9{uq — vq)G'^{u,v) . (5) 

where Gp and Gp are the chronological (Feynman) and anti-chronological (anti-Feynman) 
Green's functions, respectively, and G^ and G^ are the (quantum) Wightmann distribution 
functions. Similar decomposition can be done for the contour self-energy Sc- By using the 
2PI effective action techniques the self-energies on different branches are obtained by the 
following functional differentiation (see eg. [T2]): 



^afe/ , .X _ , , ^^2[G] 



^^\u,v)^-tab ; ' (6) 

where r2 is the sum of all two particle irreducible vacuum graphs of the theory and the 
indices a, b refer to the position of the arguments u and v, respectively, on the complex 
Keldysh time path (for example a = -|-1(— 1) implies that u belongs to the upper (lower) 
branch of the time contour in Fig. [TJ) We use the same notation S^ = — S+^, etc. as 
for propagators ([5]). Using these definitions and the relations: /^ d^n -^ Yla^ 1-00^"^^ ^"^^ 
^ci'^'^^v^) — > a6abS{u^ — v^) between complex Keldysh time and the usual real time variables, 
we can write Eq. (j4]) in the following matrix form: 



G-i ® G = ^3 5 + S ® a^G, (7) 



where 

G 



Gp -G< \ ^ / Sf -S< 



rp — Ur \ V — I ^ 

and (T3 is the usual Pauli matrix, and we defined a shorthand notation ® for the convolution 
integral: 

f®g = j d^zf{u,z)g{z,v). (9) 

We have also left out the labels u and v where obvious; for example 5 = (5^(u — v). 



2.1 KadanofF-Baym equations 

It's appropriate to further define the retarded and advanced propagators (a similar decom- 
position obviously holds for the self energy function S): 

G"^(u,tO = -0(t,O-n°)(G'< + G'>). (10) 

Equations ([7]) take on a particularly compact from when written in terms of these new 
Green's functions: 

{Gq^ - Z''^) ® G''"" = 6 (11) 

(G^^ - SO ® G'<'> = S<'>®G". (12) 

Equations ( TTTl) and ( IT2|) are called pole equations and Kadanoff-Baym (KB) equations, re- 
spectively. In general, the former will fix the spectral properties of the theory, while the 
latter will give the dynamical evolution, including quantum transport effects. Indeed, in 
the classical limit the KB-equations flT^ will reduce to well known Boltzmann transport 
equation for the phase space number density p ^fTTl [T2] . 

It can be easily shown that the defined 2-point functions have the following hermiticity 
properties: 



'iG<'>{u,v)^''Y = «G'<'>(t;,u)7° 




"2G'"(m,^;)7°]^ = -iG^iv^u)^^. 


(13) 


These identities suggest to decompose G'''" into Hermitian and anti-Hermitian parts: 




Gh ^ ^(G" + GO 




A ^ 1(G^-G0 = ^(G< + G>). 


(14) 



The anti-Hermitian part A is called the spectral function. Based on (jTOi) it is easy to show 
that Gh and A obey the spectral relation: 

Gh{u, v) = -isgn(M° - v'^)A{u, v). (15) 

Since the self-energies S satisfy identities similar to flT^ . it's appropriate to define the 
Hermitian and anti-Hermitian parts of S*"'" as well: 

sh = li^^ + ^n 

r ^ l(s"-so = ^(s< + s>). (16) 



By using the definitions (1141) and (1161) it is now straightforward to show that the pole equa- 
tions ( ITTl) can be written in the form: 

{Go^-j:H)®GH + r^A = 6. (17) 

while the KB-equation (1121) for G^ gives: 

{Go' - Sj^) ® G< - S< ® Gh = ^ (S> ® G< - S< ® G>) (18) 

Assuming we will solve the spectral function from the pole equations (IT7|) . we don't need to 
consider the other KB-equation for G'^ , since from the definition (I14j) it immediately follows 
that 

G> = -G<-2iA. (19) 

Further, we know that the canonical equal time anticommutation relation of the field oper- 
ators 

{ip{t, u),ip\t, v)} = -iS^{u - v) (20) 

must be satisfied by all physical field configurations. Using definitions i^ and flT^ it is easy 
to see that relations ( l20l) imply the condition 

2A{t,u;t,v)-f'^ = 5^{u-v) (21) 

on the spectral function. This is the direct space version of the spectral sum-rule. It follows 
also directly from Eqs. ( IT5|) and ( TT71) without a reference to the commutation relation ( 120|) . 
Note that apart from a possible implicit dependence of T on G^'^ the pole equations (IT7I) are 
entirely independent of dynamical evolution. This is exactly as it should be; the dynamical 
evolution can affect the spectral solutions related to the kinematic phase space only indirectly 
by changing the ambient conditions in the plasma in which the particles are moving. 

2.2 Lagrangian density and the mixed representation 

We now want to write the spectral and dynamical equations in detail in the mixed represen- 
tation, by Fourier transforming the 2-point functions with respect to the relative coordinate 
r = u — V. This representation is useful in dealing with the evolution equations because it 
performs a separation of the internal and external distance scales of the problem, and easily 
allows expansions in the gradients in a (slowly varying) external coordinate. At this point 
we also wish to specify explicitly the free particle Green's function Gq] in this paper we will 
consider the following CP-violating fermionic Lagrangian 

C = i^p p^ + {pLmipR + ipnm^tpL + C-mt, (22) 



where m{x) = m/j(x) +imj{x) is complex, possibly spacetime dependent mass and Ci^t is 
the interaction part to be defined later. From Eq. (l22l) it's easy to see that 



G^\u, v) = 6\u - v){i p^ -m*{v)PL - m{v)PR) , (23) 

where Pl,r = |(1 T 7^)- Next, we define the Wigner transformation of an arbitrary 2-point 
function as follows: 

F{k,x)= I d^re''''F{x + r/2,x-r/2). (24) 

where x is the average coordinate, and k is the internal momentum variable conjugate to 
relative coordinate r. Performing the Wigner transformation to Eqs. (fT7|) and (TT8|) we get 
the pole- and KB-equations in the mixed representation: 

^^+Lp^_fno-tm,^')A-e-''>{j:H}{A}-e-'0{T}{GH} = (25) 

(|^+l^^_^o_,^5^5)G'^_e-^{SH}{G^} + e-'0{r}M} = 1 (26) 

and 

iH^P, -^0 - im,^')G< - e-^^{S^}{G'<} - e"^0{S<}{G^} = Ccou, (27) 

where the collision term is given by 

Cecil = ^e-* ({S>}{G<} - {S<}{G>}) . (28) 

The <0-operator is the following generalization of the Poisson brackets: 

0{/}{^} = I [dxf ■ dug - d,f ■ dxg] (29) 

and the mass operators rho and m^ are defined as: 

mo,5F{k, x) = mRj{x)e-i^^-^^F{k, x) . (30) 

Transforming Eq. ( 12T1) in the same way gives the well known momentum space representation 
of the spectral sum-rule: 

/"^^(fc,x)7° = l. (31) 

Equations ( l25ll27l) together with the sum-rule (1311) and the identity (TT9|) form a complete set 
of equations for solving G^ and the pole functions A and Gh exactly, when the interactions 
(a scheme to compute S) and the mass profiles are specified. In practice these equations 
are too hard to be solved in their full generality, and sevaral approximations are needed 
to find a solvable set of equations. It is not clear a priori that a tractable approximation 
scheme that is general enough to treat information on quantum coherence simultaneously 
with interactions can be found. The novelty of this work is to show that such a scheme 
indeed does exist. We shall now proceed to build this scheme by first constructing the full 
phase space structure of the free 2-point functions. 



3 Free fields 

In the noninteracting case the equation (l27j) for the dynamical function G^ decouples from 
the pole functions. In this case, the Hermitian Wightmann function, defined as 

G<{u,v) = iG<{u,v)-f°. (32) 

obeys the free Kadanoff-Baym equation in the mixed representation: 

(fco +'-dt-d-{k- ^V) - 7°mo - z7V^5)g'<(A;,x) = , (33) 

which is obtained from Eq. (p7j) by setting E"* = and multiplying from both sides by 
7*^. In a spatially homogenous case the spatial gradient terms vanish and, correspondingly, 
the helicity is a good quantum number. Mathematically this follows from the fact that the 
helicity operator h = k ■ S = k ■ 7^77^, where k = k/\k\, commutes with the differential 
operator of Eq. 0331) in the homogenous limit. This fact is particularity transparent in Weyl 
basis where the gamma-matrices are given by the following direct product expressions: 

7O = pi ® 1 , a = -p^ ® a and 7^ = -p^ ® 1 . (34) 

Here both p* and a* are the usual Pauli matrices such that the p-matrices refer to the chiral- 
and cr-matrices to the spin-degrees of freedom. In this representation the helicity operator 
is just h = l®k ■ a and it's commutativity with a ■ fc, 7° and 7^ is evident. As a result one 
can introduce a decomposition of (5^ in the helicity basis: 

G<^g<®\{l + hk.a), (35) 

where g"^ are unknown Hermitian 2x2 matrices (for /i = ±1) in chiral indices. When this 
decomposition is inserted into Eq. (l33l) one obtains an equation 



(k^+'-dt)g^ = Hg<, with H = -h\k\p' + m^p'-m^p^. (36) 

Even this equation would be impossible to solve exactly, because the mass operators mo,5 
involve gradients to arbitrary orders. In the mean field limit the gradients drop out however, 
and H becomes a local Hermitian Hamilton operator 



H ^-h\k\p'' + mRp'+mip'=\ II ^ -., \=H. (37) 




One can always decompose the equation ( 1361) into two distinct sets of equations based on 
hermiticity. In the mean field limit, where H = H\ the Herm^itian (H) and anti-Herm,itian 
(AH) parts become: 

(H): 2kog< = {H,g<}, (AH): tdtg< = [H, g<] . (38) 



respectively. Note that (H)-equation is an algebraic matrix equation because it does not 
contain any time derivatives. It will give rise to phase space constraints for the components 
of (7^, as will be shown in detail in the following subsection. The (AH)-equation in flH^ 
on the other hand contains an explicit time derivative of g'f^ and is thus called "kinetic 
equation". This equation clearly has the standard form of the equation of motion of a 
density matrix in the Schrodinger picture. Indeed, it is easy to see that in the homogenous 
case the Dirac equation for the wave function ip = {Lh, Rh) (displaying only the nontrivial 
chiral components) of a given helicity h becomes just 

idt^J = Hi; , (39) 

where H is given by the mean field limit Eq. flH7|) . Given Eq. flHUjl the equation of motion 
of the form (AH)-equation in ( l38l) follows immediately for p = ipip'^ . The Hamiltonian H 
clearly has eigenvalues corresponding to free particle and antiparticle states: /cq = icu/c = 
±(P + |?T2p)^/^, and Eq. ( l39l) describes the mixing of these states in the case of a time 
dependent mass term. Despite the apparent similarities to the Dirac equation approach, 
equations (!38|) are mathematically very different from Eq. (139|) . Indeed, it turns out that 
the (H)-equations impose a singular shell structure for g^ that prevents us from interpreting 
it directly as a physical density matrix and which consequently makes the (AH)-equation 
meaningless before a sensible integration procedure is defined. 

3.1 Shell structure 

The (H)-equation ( l38l) is most conveniently discussed in the Bloch-representation for g^: 

9t ^ \ {g'o + alp') , (40) 

where g^ are real- valued functions, because of the hermiticity oig'^. In the Bloch-representation 
( HOj) the (H)-equation (l38|) decomposes into the following four real-valued "constraint equa- 
tions" (CE): 

kog^ + h\k\g^-mRg'l + mig^ = 

kog^ + h\k\g^ = 

kogi-mngl; = 

kog^ + mjg^ = 0. (41) 

This set of linear homogeneos equations can be written as an equation B^pg^ = 0, where 
the coefficient matrix is (index ordering is here defined as a = 0, 3, 1, 2): 

/ ko h\k\ —mji mj \ 



B 



h\k\ fco 

—rriR /co 

\ mi fco / 



9 



(42) 



A homogeneous matrix equation may have a nonzero solution only when the determinant of 
the matrix vanishes. Here the determinant is easily evaluated to give: 

det{B) = (kl -P- \m\^) kl , (43) 



The vanishing of the determinant (1431) yields two different classes of spectral solutions with 
different dispersion relations. First, there are the usual mass-shell solutions with k^ — k^ — 
|mp = 0, but we find also new /cq = 0-shell solutions. These solutions turn out to be the way 
the quantum coherence effects are introduced in the present approach. Let us next examine 
the matrix-stucture of these spectral solutions. 

3.1.1 A;o 7^ -solutions; free particle mass-shell 

Let us first assume that k^ ^ 0. Then the constraint equations (1411) clearly have the solution: 

93 = -^9o^ ^i=-7-^o, 92 = —rr9o^ (44) 

Kq Kq Kq 

and 

{kl - P - ImH^o' = 0. (45) 

Equation (j^3|) has the spectral solution 

g^,{koAk\;t) = ATTf^ {\klt)\ko\d{k'-\m\') 

= 27rf^^^{\klt)6{ko-Sk,^k' + \m\'), (46) 

where Skg = sgn^ko). So these solutions indeed live on the energy- momentum mass-shell 
corresponding to the dispersion relation 



ko = ±uJk = ±^Jk^ + \m\^. (47) 

Now, using (I44j) and (I46|) we can write the mass-shell contribution for the full chiral g^-matiix 
as follows: 

This solution clearly describes either a particle or an antiparticle eigenstate of helicity h and 
momentum k. 
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3.1.2 ko = 0-solutions 



Setting /co = in the first place we find out that equations (HTll have a new class of solutions, 
which obey the relations 



9's 



\k\ \k\ 



9'o 



0, (49) 

while the components gi2 are unconstrained. The corresponding spectral solution is then 



ghfi-sikoAk\;t) = 71 



fiilklt) 



+ f2m,t) 



hmii/\k\ 
1 



1 

-hmn/\k\ 



-hmj/\k\ 



hmi/\k\ 



S{k 



0), 



(50) 



— * — * ^^ 

where fi{\k\,t) and f2{\k\,t) are new unknown real functions. These solutions live on the 
ko = 0-shell, and they cannot be related to particles and antiparticles alone, since those 
should have energies ko = ±ujk respectively. On the other hand one would expect that the 
density matrix should somehow contain the information of the quantum coherence between 
the particles and antiparticles, just as does the corresponding Dirac equation. Thus we make 
a natural interpretation: the additional ko = 0-solutions [5U) describe the quantum coherence 
between particles and antiparticles with same helicity h and opposite momenta. 

Combining the solutions ( HHl) and (l50l) gives the most complete solution that satisfies the 
constraint equations (HTj) for a given helicity h and momentum \k\: 



g^{ko,\k\;t) 



9h,m~sikoAk\;t) + gfo^^{ko,\k\;t) 



27rf^^^i\klt)\ko 
fi{\k\,t) 



— h\k\/ko m/ko 
m*/ko l + h\k\/ko 



Sik' 



\7n\ 



vr 



hmii/lkl 
1 



-hmR/\k\ 



f2i\klt) 



-hmj/\k\ —i 

i hmj/\k\ 



S{k 



oj, 



(51) 



This solution contains both particle and antiparticle states as well as their coherence in 
separate singular shells in the phase space, as we promised. The solution fl^ is a distribu- 
tion which is best understood as a functional phase space measure, parametrized by weight 
functions that must be the real physical objects whose evolution we are interested in. It is 
clear that setting the form fl^Tl) directly into the (AH)-equation in flH^ does not lead into a 
sensible equation of motion however. Before discussing this problem further, we shall first 
discuss the spectral structure of the pole functions A and Gh- 



11 



3.2 Spectral function and Gh 

In the noninteracting case the equation of motion for the spectral function A Eq. (l25l) is 
identical to that for G^. As a result, the most general solution for A is of the form of 
Eq. dSSD: 

A-i"" = Y^ah®]-{l + hk-a), (52) 

h 

where the chiral matrix a^ is identical to the most general solution fIFT]) for g^, with four yet 
undefined spectral on-shell functions /^-^ for both helicities. However, the spectral function 
must in addition obey the sum-rule Eq. (IHT]) . This is enough to completely fix the values of 
all four on-shell functions: 

fi-^=\ and /l^/ = 0. (53) 

for both helicities. With these values the full solution for A becomes: 

A = 7rsgn(A;o)(^ + mji — i'j^mi)6{k'^ — \m\'^) . (54) 

This is just the familiar result for the spectral function in thermal quasiparticle limit (see for 
example [12] )• The spectral function is now completely determined, and it doesn't contain 
any dynamics at all. Especially, it does not have any contribution from the ko = 0-shell 
describing the coherence between particles and antiparticles. This is what we should expect, 
since coherence is a dynamic phenomenon, so it should not show up in the measure of the 
one-particle phase space. Moreover, it should be vanishing in the statistical equilibrium 
limit. 

For completeness we consider also the pole-function Gh in the noninteracting case, al- 
though we shall not need this function later in the paper. From fl26l) we have 

(^ko +'-dt-a-{k- ^ V) - 7°mo - t^'^-f'm^) G^ = 1 , (55) 

where again Gh = GhIo- The only difference between this equation and the corresponding 
equations for G^ and ^7° is the factor of unity on the right hand side. So we can again 
make the decomposition 

GH = Y.g'ii®\{'^ + hk-a), (56) 

h 

rewrite ( l55l) in component form and separate the equations to constraints and kinetic equa- 
tions. The only difference in these equations for g^ with respect to those for g^ occurs in 
the first constraint equation, which now becomes 

kogm + h\k\9m - mR9Hi + ^i9h2 = 1 • (57) 
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All the other component equations are identical to those for g^ in (HTj) . It is straightforward 
to solve g'^i in terms of g'^Q, and putting these back to ( 1571) gives 

{kl-k'-\m\')gj,, = l. (58) 

Now, if we were careful in keeping the e-factors through our computation, we would have 
found that the function g^^ solving the equation (158|) is in fact the principal part distribution 

9m = PP ,2 \ 12 • (59) 

It is again straightforward to show that the corresponding 4 x 4-function is given by 

G^ = (^ + m^ - z7^m,)PP ,, ^. „ . (60) 

K, \Tfl 

This form saturates the (momentum space equivalent of the) spectral relation ( IT5l) between 
Gh and A and hence is the complete solution for Gh in the noninteracting mean field limit. 
The solutions (13^ and (lUUI) guarantee that also the retarded and advanced propagators C''^ 
do not contain any contribution from the coherence solutions. 

3.3 Equilibrium limit for G^'^ 

The coherence solutions are also excluded from the dynamical functions in the thermal 
limit. First note that for free fields the equation of motion, and hence the solution, for G^ 
is identical in form to that for G^. However, the a priori independent distributions fg^'^ 
in functions G^ and G^ are constrained by the matrix relation (IT^ : G^ + G^ = —2iA. 
(Note that we drop the <, >-indices on /^-functions everywhere where there is no danger of 
confusion.) This implies that 

fsC + fs[^ = ^ a^d ft^ + ft^ = 0. (61) 

These relations hold generically as long as the spectral solutions are valid. Furthermore, 
in the thermal equilibrium limit functions G^ and G^ are related by the Kubo-Martin- 
Schwinger (KMS) boundary condition |17)i : 

G>^{t) ^ G<^{t + zf3) ^ G>^{ko) = e^''"G<^{ko) . (62) 

This matrix condition is strong enough to impose the vanishing of the coherence functions, 

fit" = (63) 



^Our sign-convention for G^ is opposite to the usual one, see Eqs. (|5l8p . hence the KMS-condition does 
not involve an explicit - sign. 
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and setting the mass-shell distributions to 



f^< = n^^iko) 



/",> = l-neq(/co), (64) 

where rieq(^o) = l/(e'^'^" + 1) is the usual Fermi-Dirac distribution. Using solutions (16311641) 
in Eq. (15T1) . and summing over the helicities in the decomposition given by Eq. (l35l) one finds 
that: 



iG^^ = 27rsgn(/co) {jt + m^ - iYmj) ncq{ko) S{k 



ml 



iG> = 27rsgn{ko){lt + mR-i-f^mj){l-n^^{ko))6{k^-\m\^), (65) 



which are recognized as the standard thermal equilibrium propagators [T7]. Generally, in 
our treatment, the functions /^^(|A;|,t) are time-dependent and carry information of both 
quantum coherence and of statistical out-of-equilibrium conditions. 

4 Weighted density matrix 

The dynamical evolution of a free system should be described by the kinetic (AH)-equation in 
fl38|) . but the singular structure oi g^{kQ, \k\;t) complicates the matters. Since distributions 
are only well defined inside an integral, the equation for g'f^ must be integrated one way or 
the other to get sensible evolution equations for the on-shell functions /±,i,2; which are the 
objects that must carry the physical information about the system. The necessity of such 
an integration is actually something to be expected, because we can never have a complete 
information about the variables that define a system under consideration. To quantify this 
thinking, we introduce the physical density matrix as a weighted average of the original 
distribution matrix g^ [7]: 

p^^iko,\k\,h;t) = J2j -^j^, Wiko,\k\,h \\ k'„\k'\,h' ; t) g<{k'„\k'\;t) , (66) 

where the weight function W{ko,\k\,h \\ kQ,\k\',h'; t) encodes our knowledge about the 
energy, momentum and helicity variables of the state. Our task is now to find out the 
equations of motion for this weighted density matrix, and consequently for the on-shell 
functions /±,i,2- For example, a complete information of the energy the momentum and the 
helicity of the state immediately renders the evolution of appropriate density matrix trivial 
one: dtp = 0, whose solutions are just the constant freely propagating helicity eigenstates [7j. 
For a more interesting example, where the quantum coherence effects become important, 
consider a case where we have a precise information of the helicity and size of momentum, but 
know nothing about the energy. Effectively this means that we cannot differentiate between 
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particles and antiparticles with the same hehcity and energy. This is the appropriate situation 
for computing the particle production during preheating for example, to be discussed in detail 
in the section mi below. The weight function describing this situation is just 

W^=^-J^5{\k\-\k'\)5h,H'. (67) 

Auk' 

In this case the weighted density matrix pw^: 

= j '^ g^ikoAklt) 

^ p,(|^|;t) = ((?<), (68) 

obeys the evolution equation of same form as g^: 

dtPh = -i[H,ph]. (69) 

Here the commutator [H, ph\ does not vanish, giving rise to a nontrivial time dependence for 
the nonsingular weighted density matrix ph- In particular, the ko = 0-shell functions /i,2 are 
now directly related to the components of the density matrix and affect the evolution of the 
mass-shell functions f± as well. Using the Bloch-representation ph = |((fi'o) + io^) ■ <?)) we 
have the following relations between (qq) and /±,i,2: 



(^o) = 


= /H/' 


(^l) = 




(^2^) = 




(^3) = 


UJ 



^A'--^m, (70) 

\k\ \k\^ 

where u = a/Zc^ + |mp. Note that the with the weight fl67|) . the physical density matrix p^ 
is just the zeroth moment of the original distribution matrix gf^ with respect to the energy. 

5 Physical quantities 

Equations (l7Up can be inverted to obtain f!^i2 in terms of moment Bloch functions {g'^)- 
We are especially interested in the expressions for the mass-shell distributions f±, defined 
in Eq. (HUj) . which in our approach are directly related to the particle and antiparticle 
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number densities. Indeed, according to Feynman-Stuckelberg interpretation the phase-space 
particle number density is n = /+, while for fermionic antiparticles n = 1 — /_. With these 
identifications, using the inverse relations of Eq. (170]) we find that for a given 3-momentum 
k and helicity h the out-of-equilibrium particle and antiparticle numbers can be written as: 

Hh = ^{-h\k\{g',)+mM)~mj{g!^))-^{g!;) + l. (71) 

Setting Tr(p'') = (g^) = lo these expressions reduce to the ones obtained in ref. [2Jj, where 
they were derived using the solutions to a Dirac equation and a Bogolyubov transformation 
to diagonalize the fermionic Hamiltonian. Indeed, this definition of the particle number 
in terms of the moment functions was one of the main results of the paper [2]. Here the 
definition of the particle number is trivial, and we introduced the expressions (1711) merely to 
show that our definition does agree with the other approach. In a future work [8], we shall 
show that in the bosonic case our particle number differs from the one obtained in ref. [2], 
but is consistent with the definition of the ref. [IH] . 

It is also interesting to see what kind of expressions other physical quantities like energy 
density and pressure will have in terms of the components (f?^) or the on-shell functions 
f±,i,2- These quantities are defined as the ensemble averages the energy momentum tensoio 

e^"" = J Ui-i^'d"^ - d^^Y'^^ +fi^u (72) 

and so, for example for the energy density we get 

d^k 

d^k 



{n{x)) = {e'\x)) = I j-^^Tr[{i.k + mn + tmn'yG<ik,t) 

h\k\{g'3) + mH{g1)-mr{g',) 



h 

E 



(27r)3 

d^k 

(2^ 



■ u 



fcKh+^fc/^-l) • (73) 



According to expectations the result is simply a sum of free particle and antiparticle contri- 
butions. The last term in the last row is the sum of corresponding vacuum energies. For the 
pressure we get instead 

(P(x)) = (r(x)) = |^Tr[yF.G<(fc,t)] 



■^Physically this constraint corresponds to setting the chemical potential to zero. Indeed, since (g^) = 
f+ + f- = n^i^ — njr^j + 1, we see that setting (gg) = 1 reduces to ng^ = n^j^. 

'^Different signs of terms involving h and mj are due to a different convention in our definition of the 
Hermitian Wightmann function Eq. (j32p . 

^Here we use the symmetric (Belinfante) version of energy-momentum tensor |19j . However, the same re- 
sults would have been obtained using the canonical tensor T^^'^ = ^ [ij'^d'^ — g'^" (ij'^d^ — mu — 17^7717)] ^. 
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Now we see that in addition to normal free particle and antiparticle terms there is an explicit 
contribution from the coherence shell functions fi^2, signalling that at the quantum level 
the pressure is different from the statistical one. In most cases the quantum effect would 
be unobservable however, since the coherence functions /i^2 are typically oscillatory with 
microscopic time-scales Atosc ~ l/^i so that the classical thermodynamical pressure arises 
from Eq. fl74p when it is averaged out over any time-scales exceeding the quantum scale 

6 Interacting fields 

Having set up the density matrix formalism for treating quantum coherence phenomena in 
classical backgrounds, we now wish to extend our work to include interactions. To this end 
we must use the full Kadanoff-Baym equations fl25ti27p . In their complete generality these 
equations couple nonlinearily the three Green's functions G^, A and Gh- Solving these 
equations simultaneously would be an overwhelmingly difficult problem however, and we 
shall adopt a series of approximations to extract the relevant physics in what becomes the 
quasiparticle limit. The key idea in our approach is to divide the problem into two parts. 
We first need to find a reasonable approximation for the phase space of the problem in terms 
of on-shell distributions as we did earlier in the case of the free fields. Second, we must 
find the equations of motion for the on-shell distribution functions including interactions. 
Looking at equations ( l25ll271) one immediately sees that the couplings between equations due 
to terms involving the pole function Gh are causing most problems along the way to any 
Boltzmann-equation type approximations for the problem. If these can be neglected, the 
equations for G^ and A will decouple and one does not need to solve Gh at all. Fortunately, 
as we shall see, this indeed is a reasonable approximation in the weak coupling limit. As a 
result of this approximation A and G*^ will continue to have on-shell solutions, which allows 
us to go through the procedure leading to our density matrix formalism, yet including the 
effects of decohering interactions. 
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6.1 Quasiparticle approximation 

Let us first consider the weak coupling approximation for the pole equations ( I25ll26ll . To 
zeroth order in gradients, but for arbitrary Ti^ and F we obtain: 



1 + (G„r)2 

where Go now includes the real part of the self-energy S/^: 

Gq^ = jt-iTiR- imn^ - T,H , (76) 

and r is defined in equation ( IT6l) . We see that A no more has a spectral solution, and the 
phase space is truly 4-dimensional. However, A clearly reduces to a spectral form in the 
limit r ^ C@: 

A -^ 7rsgn(/co)<5(G^^) 

= 7rsgn(fco)Godet(Go-^)<5(det(Go-i)), (77) 

where the determinant inside the delta-function now gives rise to a modified dispersion 
relation: 

det(|^ -niR- imn^ - Eh) = 0. (78) 

For the pole-function Gh we would obviously get the corresponding principal value as in 
noninteracting case (see section [3l2!) . The limit of taking F — *> 0, while keeping Eh finite is 
just the well known quasiparticle approximation. It is often a reasonable approximation to 
be made in the weak coupling limit. Technically one should require that Eh is of lower order 
in the coupling constants than is F, and this often indeed is the case: for gauge interactions 
for example one finds Eh ~ g^ and F ~ (7^ in the lowest order in the gauge coupling g. 
However, even when the coupling hierarchy is not there, the quasiparticle limit can be a 
useful first approximation. 

6.2 KB-equation with collisions 

We now turn our attention to the dynamical equation ( 1271) . which in the mean field limit 
becomes: 

(^ _ 1 ^^ -m^ - tmn' - Eh) G< - E<Gh = Ccoii , (79) 



^Note that F is in fact a 4 x 4 matrix operator, so with the hmit F ^ we mean that the couphng y ^ 0. 
Actually wc have to keep a finite, but arbitrarily small imaginary part in the definitions of propagators, due to 
the inclination of the Keldysh path in the complex time plane. Thus in the collisionless limit F — > sgn(fco)£l4- 



where 

Ceoii = i(S>G<-S<G'>). (80) 

Clearly the Gjif-niixing term prevents Eq. ( 1791) from providing unique solution for G^ even 
when the spectral function is known. However, neglecting this term is in fact consistent 
with the quasiparticle approximation. One can see this from the fact that Eq. ( l25l) for the 
spectral function A can be obtained as a sum of the evolution equations for G^ and G^ , 
and that the mixing term ~ TGh in this equation arises from the sum of G/^S^'-^-terms 
in the equations for G^'^ . Since neglecting the F- mixing was precisely what defined the 
quasiparticle approximation for the pole-equations, corresponding terms should be discarded 
in the equations for G^'^ as well. Thus, in the quasiparticle and mean field limit we have 

{jt-'-p., -niR - zmn' - ^h) G< = Coll . (81) 

The collision term can be written in terms of F and A as follows: 

Ccou = -«FG< + iJ:<A , (82) 

so that given a solution for A equation (IHT!) can be solved in weak coupling limit to give 
G^. The overall strategy is the same as before: one first divides equations into kinetic equa- 
tions and constraints. Consistency with the quasiparticle approximation requires that one 
must neglect all collision terms in the constraint equations, since including them would give 
corrections that are of the same order as terms neglected in quasiparticle approximation. 
When this procedure is followed the constraint equations for G^ lead to the same quasipar- 
ticle mass-shell solutions for G^ as for the spectral function, plus the additional solutions 
describing the coherence at /cq = 0-shell, because for G^ the latter are not suppressed by the 
spectral sum-rule. However, our main interest in this paper is to study the effects of decoher- 
ing collisions on the quantum coherence of the system. Because these do not qualitatively 
depend on the modifications to dispersion relations, we will set Hh to zero for simplicity in 
what follows. This obviously reduces the functions A and Gh to their coUisionless limits 
given in Eqs. (I54|60p . and the kinetic equation eventually becomes simply 



(A;o +'-dt~a-{k- ^ V) - 7°mo - 27°7'^5)g'< = 27°Ccoii7° • (83) 

Until now our analysis has been entirely independent of the particular type of interactions. 
However, it is not possible to analyse the problem further before explicitly defining the 
interaction terms. Before doing that we shall notice that whatever the form of the interaction, 
if the self-energies S^'^ are thermal obeying the Kubo-Martin-Schwinger (KMS) relation 
T,^{k,x) = e'^''°T,^{k,x), the collision term reduces to: 

Ccoii = -^rKMs(G<-G<) . (84) 
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This follows from the KMS-relation combined with the relation iGf = 2neq(/co)^ given 
by Eqs. ( 15^ and (l65l) . The form of collision term ( 18^ is familiar from relaxation-time 
approximation and it yields the relaxation to thermal equilibrium in time scale l/F in the 
absence of any driving terms. 

After the interactions are specified and an explicit form of F is known we can proceed 
to write down our equations of motion. These equations will be generalizations of the 
equations ( l38l) in the free field case. As we mentioned above, the constraint equation (H) 
will be untouched because we are working in the quasiparticle mean field limit and we 
neglect the term S// for simplicity. The kinetic equation (AH) will receive contributions 
from the collision terms however. It is very important to realize, and this can already be 
seen from equation ( l83l) that as the singular solution for the function G^ is introduced into 
the kinetic equations, we will encounter projections of the collision term F on all different 
shells contributing to G^. That is, the external momentum configuration entering to the 
evaluation of F will depend on the particular shell multiplying it in the collision integral. In 
particular, the coherence shells will pick interaction terms with T{kQ = 0), which in general 
are completely different from the terms involving the usual mass-shell functions. 

6.3 Computation of the self-energies 

In this subsection we give an explicit evaluation of the self-energy corrections arising from 
a specific form of the interaction. This is a rather technical calculation whose results are 
particular to the chosen interaction. A reader not interested in these details can skip this 
subsection and move directly to the next section where we will continue developing the 
kinetic equations given the results from this subsection. We choose to consider the following 
left-chiral non-diagonal Yukawa interaction term: 

>Ci„t = -y V^L^ Qr + h.c. (85) 

where '0 is the considered fermion (quark) field, q is some other fermion field and </> is a 
complex scalar field. As mentioned before we use the two-particle irreducible (2PI) action 
method to calculate the self-energies (E]). The lowest order 2PI-graph based on interaction 
f lS^ is presented in Fig. EK, and it gives the contribution 

Tspi = -\y\^ I d% d% Tr [PrG,{u, v)PlG{v, u)] A(n, v) , (86) 

Jc 

where G, Gq and A are propagators of the considered fermion, quark and scalar field and 
the integration is over the Keldysh path. From this we get fermion ip self-energies: 

E'^\u, v) = -rah^^^ = i\y\'PnGf{u, v)PlA'^\u, v) (87) 
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Figure 3: Diagrams contributing to the 2PI-action and the self-energy at the one loop level 
due to interaction fl85l). 



and in particulaiu 



T.<>{u,v) = i\y\'PRG<^>{u,v)PL^<^>{u,v). (89) 



We obviously need these quantities in the mixed representation and after a Wigner trans- 
formation the self-energies become: 



E<'>(fc, x) = i\y\^ j ^^P^G<'>{k', x)PLA<'>{k - k', x) . (90) 

To keep things simple, we will now assume that the quark- and scalar distributions appearing 
in the loop are thermal. The appropriate quark propagators (with real constant mass) can 
then be read directly from (l65ll : 



^Gl,^{k) = 27i{jk + mg)6{k^-ml)sgn{ko)n,^{ko) 

iGl^^{k) = 27r(^ + m,)5(P-mJ)sgn(A:o)(l-neq(M), (91) 



while the equivalent expressions for the bosonic propagators are [12] 

zA<(A;) = 27:5{e - ml)sg^{ko)nt^{k,) 

tA>^{k) = 27r5(A;2-m5)sgn(fco)(l+<(M), (92) 

where the thermal fermion and boson distribution functions ricq and nf are 

^-^(^o) = ^^^^r^ ' ^t(^o) = ^^,, _ ^ • (93) 

In the thermal approximation the explicit x-dependence of the propagators G^'^ and A^'^ 
drops, and it might appear that also S^''* then becomes x- independent. That this is not so 
follows from the nontrivial dependence of Af'-^{k — k') on the exteral 4-momentum k. When 
k enters to the on-shell delta-functions it introduces dependence on the mass of the external 



^The 2PI-formalisni is not necessary for obtaining these results. Equivalently with Eq. (|86|) one can 
directly write down the self-energy with complex time variables: 

^c{u,v) = i\yfPRGliu,v)PLAc{u,v), (88) 

from which (|87p readily follows by making the appropriate choices for the position of variables on the time 
path. 
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field ijj. Computing S^'^ is now a simple matter of substituting tliermal propagators 
to tfie expression for the self energy ( l90l) and using the delta-functions to perform as many 
of the integrals as possible. Before doing the actual computation let us note however, that in 
the thermal limit S^'^ can only depend on two independent 4-vectors: k^ and the 4-velocity 
of the plasma m^. Given the chiral structure of the interaction the most general form for any 
E"^ is thus 

E"''(fc) = ( A"^ ^ +5"^ i) Pl, (94) 

where the isotropy of the thermal distribution implies that functions A""^ and B""^ can only 
depend on ko, \k\. Note that at this point S*^^ is actually independent of the external 
coordinate x; the x-dependence (actually only t-dependence in this paper) is only introduced 
through projections to mass shells where the dispersion relation ko = ±ujk{x) depends on x 
through the mass of the ip-field. Note also that Eq. (!M|) is more general than the specific 
interaction we are studying here: it is valid for any L-chiral interaction and to any order in 
the loop expansion as long as the loop particles are assumed to be thermal. Substituting the 
expressions ( l9Tll92l) to Eq. ( l90l) we get, after integration over k'^: 

d^k' 

ne,i±u^l)ni{koTu^l)^'^PL 



ij:<ik) 



2My\'Y.I 






UJ 



k-k> 



(95) 



where a;^ 



•p^ + "^q<A ^^^ ^± — (^fc'' i^O- The angular integrals in the expression ( l95l) 
can be further evaluated in spherical coordinates with the result: 

S<(A;)= [S< 7° - S3< (L 7)1 Pl , (96) 



where k = k/\k\ and iS^3 = i^Q^^ko, \k\) are real-valued functions of the phase space. 
For our present analysis we need to evaluate the self-energies S^'-* both on the mass-shell 



k^ — k"^ = |m(a;)p as well as on the k^ = 0-shell. On the mass-shell we get: 



zE^(A;o = ±uJk{x), \k\) 
iT.f{ko = ±uJk{x), \k\) 



2rr2 



\yn 

8n\k\ 
87r\k\ 



\h{ko,\k\)\ 



^[\hikoAk\)\ 
\k\ 



\a\ \7n\ 



|/o(A;o, 1^1)1 



with 



and 



^"(^-'^')-^(^)j(_/^^ ^e. + l)(eV.-.-l) 



(97) 
(98) 

(99) 



Imp + m^ 



a 



P 



ml ko ^sE^ko 



2|?Tip T |mpT 

X'l''{\m\\mlml) \k\ ^ v/ip,|fc| 
2|??7,p T |?npT 



(100) 
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where cufc = a/A;^ + |r7ip and £'=„ and p* are the energy and momentum of the decay prod- 
ucts in the decay frame, ^i is the invariant mass of the decaying (heaviest) particle and 
X{a,b,c) = (a + fe — c)^ — 46c is the usual kinematic phase space function. The peculiar 
appearance of the Heavyside step-function with the argument A = A(|mp,m^,m^) causes 
these expressions to automatically take care of the correct mass hierarchy of the fields; the 
mass m = m{x) may be varying in spacetime, so ip can for example change from being the 
lightest to the heaviest of the three fields. On the ko = 0-shell we get instead: 

iYl<{ko = 0,\k\) = ^^ / dy . ;, , (101) 

°^ ' '^ 87r|A;| J ^^ smh y ^ ^ 

I I 2|fc|T 

^S3<(A;o = 0,1^1) = 0, (102) 

Now A > always so there is no need for an explicit step-function. One crucial difference 
between the mass-shell and the /cq = 0-shell self energies is that the latter are completely 
x-independent while the former may be x-dependent, because of the possibly x-dependent 
mass, as mentioned earlier. 

Because we are computing E'^'"s in the thermal limit, the expression for S^ can be 
obtained from that for S^ by use of the Kubo-Martin-Schwinger (KMS) relation: 

S>(A;)=e'^'^»S<(A;). (103) 

This relation can be seen to emerge from (1951) . where we should obtain S^ -^ S-* by replacing 
n| — »■ 1 — raj^q and nf^ —>■ 1 +nf^. KMS-relation (I103p follows then immediately by use of the 
relations 1 — ncq{ko) = e^'^°?7,cq(A;o) and 1 -|- nf^{ko) = e^^°nf^{ko). Self-energy components 
Sq 3 of course obey the KMS-relation separately. Moreover, we find that 

T{k)='-{l + e'"'^')j:<{k). (104) 

Expressions (I96p and (11041) complete the computation of all required self-energy functions 
needed for the further analysis of Eq. fl83|l . 

6.4 Evolution equation for the density matrix 

Using the particular form of the self energy ( l96l) and the KMS-relation (I104p . which yield a 
thermal collision term of the type Eq. (1841) . and the helicity block-diagonal decomposition 
(l35l) for G^, the equation ( l83l) reduces to the following matrix equation for the chiral part 



zD{g<~{g<U), (105) 



9h- 






(ko + \dt-H)g< = -i. 


where H = 


— h\kp^ + mf>p^ — mip"^, as earlier. 




D = ^{l + Ps)Tf, with 
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r^ = To - /lEs (106) 



and {g^)cq is the thermal equihbrium hinit of g^ defined in Eq. fl5ip . Taking the Hermitian 



and anti-Hermitian parts of equation fllOSp . and neglecting collisions in the Hermitian equa- 
tion, consistently with the quasiparticle approximation, we find a generalization of equations 
( l38l) to the case with collisions: 

(H): 2kog< = {H,g<} 

(AH) : d,g< = -z[H, g<] - {D, g< - {g<U} . (107) 

The only difference to the free field case is the appearance of an anticommutator in the anti- 
Hermitian equation containing the interaction matrix D. The precise form of the matrix 
D in chiral indices will depend on the form of the interactions, but the generic form of the 
equations in 01U7p including an anticommutator with a generic matrix D is universal to the 
quasiparticle mean field limit when the self-energy is computed in thermal approximation. 
The simple form of D in Eq. 01U6p refiects the chirality of the interaction: (1 + P3)/2 is just 
the 2-dimensional version of the left- chirality projector. 

Just as in the noninteracting case, the (AH)-equation is ill-defined without an integration 
procedure. Here we are interested in the evolution of the weighted density matrix relevant for 
the particle production during preheating in the homogenous time-dependent background. 
The momentum and helicity remain to be good quantum numbers even in the presence of the 
interactions, and so we can take our physical density matrix to be diagonal in these variables. 
Furthermore, since we have no a priori information of the energy of the relevant solutions 
we impose a fiat weight on the energy variable. That is, we integrate the (AH)-equation in 
(TT07I) with the weight ((671) introduced in the section H Wi = ^^6{\k\ - \k'\)6hh'- After 
some manipulations we get the new equation of motion for the physical density matrix ph 
defined in Eq. ( l68l) including interactions: 

^tPh = -^[H,p^,]-Ig, (108) 



where 



h = 



l^{D,g<-ig<r} 



mO 



- hVr 



rn* ( f /cq\ 

m* ( f /eq\ 

2^1/0 - Jo ) 



tih- 


/D 







iS(^o - 


~ /o 



i(/l+^/2) 
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where we use shorthand notations /o = /+ + /_ and /s = f+ — f-, and similarly for /q 3 with 
f^ = neq{±uj). Yet we have denoted k = \k\, and we have dropped the helicity index h in 
the superscripts of all fa for convenience. The F-functions appearing as coefficients are now 
having contributions from different shells on the phase space as we predicted at the end of 
the section I6.2[ We have used the definitions: 

rm(o,3)(|^M) = ^0,3(^^0 = i^it),\k\) (positive mass-shell) 

Tood^l) = ro(A;o = 0, 1^1) (A;o = -shell) . (110) 

These are the only independent functions, since F on the negative mass-shell is related to 
that on the positive with Fo,3(— /i;o, |^|) = ±Fo,3(/i;o, |^|)- Further, F3 vanishes on the ko = 
-shell. It is because of the different values of F-functions in different phase space shells, 
that the integrated equation of motion (llOSp cannot be written in an equally simple matrix 
form as the original one (I107p . We have written the collision integral (I109p in terms of 
the on-shell functions fa to get the simplest possible expression. Of course, in order to 
perform any practical calculations one has to use the relations between pij and fa, which 
can be found using the relations (170]) . In fact, in our numerical calculations in sections 17.11 
and 17.21 we found it easiest to express all quantities in terms of the Bloch-components ((7^). 
Again, the form (llOSp is generic for a density matrix defined by the weight 06 7p . but the 
precise form Ig strongly depends on the interaction; the more complicated the interaction, 
the messier Ig becomes. However, there is always a well defined matrix expression for Ig 
which in the thermal limit is linear in the on-shell functions, involving projections of the 
interaction matrix components on the positive mass shell and the ko = 0-shell. 

The generic features of the thermalization and decoherence effects due to collisions are 
a bit obscure in Eq. (llOSp because of the complex matrix stucture of the collision integral 
(I109p . They become much simpler if one looks directly at the evolution equations of the 
on-shell functions fa in the constant mass limit. Using the relations (1701) to express pij in 
terms of fa, we find that in this limit: 

dtf± = -(l-/i^)F^o(/±-/r) + --- 

dtfi = -2/i[-A:/2 + m,(^/i-^/2)] -roo/i + ... 



dj, = -2h\kf\+mn{'^f^-'^f,) 



roo/2 + ..., (Ill) 



where we have written down only the diagonal interaction terms. Since all of those are 
negative (Fmo,roo > 0), we can immediately see the tendency of interactions to thermalize 
the mass-shell functions by setting f± — »■ f^ as well as to give rise to decoherence /i^2 ~^ 0. 
The latter effect in particular gives a smooth damping of the quantum coherence as a result 
of collisions over a characteristic time scale I/Fqo- 
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Figure 4: Shown is the mean field number density ra^^ of produced fermions as a function of 
r = uj^pt for negative hehcity state /i = — 1 in free case (thick black line). Also shown is (thin 



dotted blue line) a function /^ = {{fiY + {f-iYY^"^^ which describes the overall amount of 
correlation between fermions and antifermions. Effects of infiaton oscillation are modeled by 
varying mass m{t) = (10 + 15 cos(2ci;^t) + isin(2ti;<^t))|/c| where lj^ = \k\ is the frequency of 
the infiaton oscillation. At r = the fermion system is taken to be an uncorrelated vacuum 
state. 

7 Applications 

Let us next study the effects of collisions numerically in two different physical examples. 
First we consider particle production at the preheating of inflation including the effects of 
collisions on the coherent particle number creation. Second, we study the approach to the 
thermal equilibrium of a quantum system including both thermalization and decoherence 
effects due to collisions. In both of these examples the 3-momentum k and the helicity h are 
good quantum numbers so that the weight function Wi deflned in Eq. ( 1671) . and the resulting 
equations OlOSHllOp with relations (1701) are the appropriate one to use in these calculations. 

7.1 Particle production at preheating 

As our flrst application we consider the particle production at the preheating stage at the end 
of the inflation, where a time- dependent fermionic mass is generated by an oscillating inflaton 
condensate. This system is appropriately described by equations (1108111 lOp . assuming that 
interactions are modelled by our left-chiral Lagrangian Eq. fl85l) . We found it most convenient 
to write the equations entirely in terms of the moment Bloch functions (g^) in the code, using 
relations (1 70 p . and by computing the particle and antiparticle numbers in the end through 
Eqs. (17T1) . This preheating problem was considered recently in ref. j2] in the coUisionless 
limit, and to facilitate the comparison, we will adopt their model for the oscillation of the 
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Figure 5: The mean field number density n^^ of fermions as a function of r for an interacting 
field (thick black line), in same setting as in Fig. HI Thin dotted blue line again shows the 
coherence function fj! giving the correlation between fermion and antifermion fields. The 
interaction is taken to be of the form computed in the section [U75| with a thermal background 



in temperature T = \k\ and parameters y = 5,mq = 0.02|A;| and ?7i^ = 0.1|/i;|. 

inflaton condensate. A simple cosine function for inflaton gives rise to the following fermionic 
mass function: 

m(t) = nio + A cos{2u^t) + iB sin(2u;<^t), (112) 

where ttlq, A, B and u^p (the inflaton oscillation frequency) are real constants. Let us first 
consider the noninteracting case with F = 0. We assume that there is no initial chemical 
potential so that rir 



kh 



n 



^^ and ((7g) = 1. This initial condition is preserved throughout 
the free field evolution, since by Eq. (l69l) one can show that dt{go) = 0. The results of our 
calculations for helicity h = —1 are presented in Fig. HI We see that the produced particle 
number increases steadily as a function of time. This increase takes place at the resonance 
peak areas while between the peaks the particle number is essentially constant. To arrive to 
this picture with a parametrical resonance needs some fine-tuning between the parameters 
of mass oscillation and the size of momentum \k\ and helicity h; indeed for the opposite 
helicity h = +1 with the otherwise same parameters we would get a completely different 
figure without a clear resonance behavior. This evolution of the particle number agrees with 
the results in ref. [2] . We show also the evolution of the coherence by plotting a function /^ = 
V/i + /I (dotted blue line) in Fig. HI We see that the particle production is accompanied 
by a steady growth of coherence. Most of the growth takes place at the resonance peaks, 
whereas outside the peaks the coherence oscillates with a constant amplitude. Already after 
a few resonance crossings the coherence has saturated to a maximum. 

Let us next consider the case with collisions, so we have F 7^ 0. Physically our model 
interaction (l85l) corresponds to decays (and inverse decays) of the inflaton generated field ip 
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Figure 6: Shown is the number density n^^ in same setting as before with changing interaction 
strengths. Thick black hne is free field case with a coupling constant y = 0. The other lines 
are interacting with y = 1 (green line), y = 3 (blue dotted line) and y = 5 (red dashed line). 

to some lighter fermions and scalars that could be taken to be the standard model particles. 
Including transport equations also to these states (and the ones they couple to), and solving 
the whole problem with the momentum dependence, one could build an entire network of 
equations necessary for a realistic simulation of the particle production and thermalization 
at the preheating. Here our goals are more modest, and we only wish to study how the 
collisions affect the evolution of the particle number and coherence. As before we will set 
the initial chemical potential to zero so that n^^ = n^^ and (qq) = 1 in the beginning. 
However, now this condition is not preserved, since conservation of (g^) is not respected 
by the collision terms. We show the evolution of the particle number and the coherence 
for a particular parametrization of the collision term in Fig. [5l Apart from including the 
interactions we are using the same parametrization as in the noninteracting case shown in 
Fig. m The difference between the interacting and noninteracting cases is quite dramatic. We 
see that with interactions the particle number drops between the resonance peaks, only to be 
regenerated again in the next resonance crossing. The drop in the particle number obviously 
results from the decay of the unstable particles, so that a steady flux of standard model 
particles (and antiparticles) is created through the resonance production and decay of the 
ip-states. Also the growth of the coherence is damped in comparision with the noninteracting 
case, and both the particle number and coherence evolution settle into a stationary pattern 
only after a few oscillation periods. In Fig. [6] we plot the particle number evolution for 
varying interaction strenghts, again for our reference set of infiaton parameters. The orderly 
parttern of the particle number evolution over the infiaton oscillation periods remains, while 
the damping effect depends on the strength of the interactions in a predictable manner. 
As we mentioned, getting the parametric resonance needs some fine tuning of the pa- 
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Figure 7: Same as in Fig. [6], but with the periodic mass term replaced by a damped oscil- 
latory mass: m{t) = (10 + exp(— 7r)(15cos(2ci;^t) + zsin(2c(;^t)))|/c|, where 7 = 0.05. Other 
parameters are again set as in all previous figures. Thick black line shows again the free field 
case and the interacting cases have the coupling constants y = 1 (green line), y = 3 (blue 
dotted line) and y = 5 (red dashed line). 

rameters. Moreover, when the resonance condition is not met, the evolution of the system 
becomes much more complicated than the orderly behaviour shown in Figs. miUl In Fig. [7] we 
plot the particle number evolution for the case where the infiaton oscillation is exponentially 
damped, leading to a new modified mass function: 



m{t) = mo + e '^'^ (A cos(2u;^t) + iB sin(2u;^t) j . 



;ii3) 



Although we are using a rather small damping parameter, the evolution of the particle 
number is dramatically changed. The orderly oscillation is replaced by an essentially chaotic 
evolution over the resonance crossings with changing conditions. The effect of interactions 
is also more pronounced, showing a very strong quantitative and qualitative dependence of 
the evolution on the strength of the interactions. 

Our results show that making the heavy particles unstable can have dramatic effects on 
the pattern of the particle number production during the preheating stage of the inflation. 
This instability is of course a necessary requirement as these states need to decay to the 
standard model particles to eventually reheat the universe. However, our formalism is well 
suited for a detailed study of this problem, if only the appropriate interactions are specifled, 
the necessary transport equation network is written down also for the daughter particles, 
and when the inflaton dynamics is also modelled numerically along with the evolution of 
the transport equation network. Moreover, our methods can be used to model coherent 
baryogenesis [1], where collisions might have large quantitative effects. Let us flnally note 
that if one introduces a CP-violation into the heavy ip-decajs, one can use our formalism to 
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directly and accurately compute the evolution of the baryon (or lepton) number, in a possible 
decaying heavy-fermion baryogenesis scenario during the particle production at preheating. 

7.2 Decoherence and thermalization due to collisions 

As our last application we consider the thermalization and decoherence of an arbitrary, 
correlated out-of-equilibrium density matrix by collisions. We consider a homogenous system 
described by the weight function fl671) . We also take the mass to be a constant here, and 
begin with a correlated initial state at t = 0, described by functions 

n(0)=n(0) = i, /i(0)=0 and /2(0) = 1. (114) 

We then use the equations fllOStlllOl) to calculate the evolution of the density matrix. We 
assumed that the decay product species q and </> are in thermal equilibrium and the param- 
eters in the interaction term were taken as follows: coupling y = 1 and the masses m = 10, 
rrig = 0.02, m^ = 0.1 and the temperature T = 10 (all in units \k\). The results of the calcu- 
lation are shown in Fig. [8l One sees clearly how the collisions smoothly damp the amplitude 
of the oscillating coherence on-shell functions /i 2- As long as the coherence functions are 
nonzero, they affect also the mass-shell functions that show a characteristic oscillatory pat- 
tern. After the coherence solutions are damped out, the mass-shell functions f± continue to 
approach their equilibrium limit, which for the current parameters corresponds to n^^ ^ 0.27, 
shown by the thin solid line in Fig. [HI This behaviour can also be seen qualitatively from 
the Eq. 01 111) in section [HI As a result the full 2-point Wigner function G^ approaches 
the thermal limit Gf given by Eq. (l65l) . as it should. The time scale for vanishing of the 
quantum coherence is here smaller than the thermalization time scale simply because for the 
parametrization chosen for the problem, the ko = 0-shell collision rate is much larger than 
the mass-shell collision rate: Tqq/Tj^ ^ 10. 

8 Conclusions and outlook 

In this paper we have developed a quantum transport formalism, simple enough to be used in 
practical calculations, for treating the coherent evolution of fermionic system in the presence 
of decohering collisions. In our derivation we used the CTP-formalism for the quantum 
field theory in the out-of-equilibrium conditions. The key element in the derivation was 
the observation that in the mean field limit and in the quasiparticle approximation in an 
interacting theory, the 2-point function has a singular shell structure, that in addition to the 
usual mass-shell solutions contain new ko = 0-solutions, which carry the information about 
the quantum coherence between particles and antiparticles. Thus the actual 2-point function 
can be understood as a dynamical phase space weight function, or measure, with a number of 
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Figure 8: Shown is the evolution of the number density ng^ (thick sohd hne) and the coher- 
ence function /^ (wiggly dotted hne) from an intially highly correlated out-of-equilibrium 
configuration towards the thermal equilibrium. The thermal equilibrium value rigq ^ 0.27 is 
shown by the thin straight line. 

unknown on-shell coefficient functions that carry the physical information about the system. 
These can be interpreted as the out-of-equlibrium particle numbers on mass shells and as 
functions measuring quantitatively the level of quantum coherence on ko = 0-shell. We 
showed that the new coherence solutions are excluded from the spectral function by the 
spectral sum-rule, and that they are also eliminated from G^ by the KMS-condition in the 
thermal equlibrium limit. The fact that coherence solutions appear only in the dynamical 
function makes sense, because they describe correlation between mass-shell states, and no 
coherence should be present without interfering physical states (Solitary coherence functions 
can however, in some cases describe pairs of localized virtual states [7].) 

We then proceeded to show how a sensible physical density matrix can be defined by 
an integration procedure which involves specifying the amount of external information on 
the system. This information affects both the definition of the physical density matrix 
and its eventual equation of motion. Finally we derived an explicit dynamical equation of 
motion for a physical density matrix corresponding to a spatially homogenous system, where 
the 3-momentum and helicity are good quantum numbers. We also computed the collision 
integrals in the thermal limit for a model interaction Lagrangian, showing in particular how 
different shells pick up different collision terms in the integration procedure, corresponding 
to the different external momentum configurations at different shells. Finally, we applied 
our formalism to the case of coherent particle production during the preheating stage of 
the inflation including a flnite decay width to the produced heavy particles. We showed 
that including collisions can have dramatic effect on the evolution of the particle number, 
in particular when one includes a nontrivial evolution of the inflaton fleld. We stress that 
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our method can be used to a quantitative analysis of the particle production during inflation 
if one models the inflaton oscillation and evolution numerically, and includes a network of 
transport equations also for the daughter states in the decay process. Finally we pointed 
out that the method is also suitable for a quantitative analysis of coherent baryogenesis and 
a computation of a baryon- or lepton number generation during decay of the preheating 
produced heavy states. In our last example we showed explicitly how an initially highly 
correlated out-of-equilibrium density matrix relaxes to a thermal equilibrium. 

We believe that the current formalism will find applications in many different contexts, 
also outside the cosmology. In particular we are currently extending the formalism also to 
the case of relativistic bosonic fields and also to nonrelativistic systems. It is also interesting 
to see how the usual fermionic flavour oscillation pattern arises in the present context when 
the formalism is extended to the case with many mixing fermion flavours [TJ [H]. 
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